Casimir Force at a Knife's Edge 
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The Casimir force has been computed exactly for only a few simple geometries, such as infinite 
plates, cylinders, and spheres. We show that a parabolic cylinder, for which analytic solutions 
to the Helmholtz equation are available, is another case where such a calculation is possible. We 
compute the interaction energy of a parabolic cylinder and an infinite plate (both perfect mirrors), 
as a function of their separation and inclination, H and 6, and the cylinder's parabolic radius R. 
As H/R — s> 0, the proximity force approximation becomes exact. The opposite limit of R/H — >■ 
corresponds to a semi-infinite plate, where the effects of edge and inclination can be probed. 

PACS numbers: 42. 25. Fx, 03.70.+k, 12.20.-m 



Casimir's computation of the force between two paral- 
lel metallic plates [l| originally inspired much theoretical 
interest as a macroscopic manifestation of quantum fluc- 
tuations of the electromagnetic field in vacuum. Follow- 
ing its experimental confirmation in the past decade 0, 
however, it is now an important force to reckon with in 
the design of microelectromechanical systems 0] . Poten- 
tial practical applications have motivated the develop- 
ment of numerical methods to compute Casimir forces 
for objects of any shape The simplest and most com- 
monly used methods for dealing with complex shapes rely 
on pairwise summations, as in the proximity force ap- 
proximation (PFA), which limits their applicability. 

Recently we have developed a formalism [1, 0| that re- 
lates the Casimir interaction among several objects to the 
scattering of the electromagnetic field from the objects 
individually. (For additional perspectives on the scatter- 
ing formalism, see references in [6[.) This approach sim- 
plifies the problem, since scattering is a well-developed 
subject. In particular, the availability of scattering for- 
mulae for simple objects, such as spheres and cylinders, 
has enabled us to compute the Casimir force between two 
spheres Q , a sphere and a plate Q , multiple cylinders [1] , 
etc. In this work we show that parabolic cylinders provide 
another example where the scattering amplitudes can be 
computed exactly. We then use the exact results for scat- 
tering from perfect mirrors to compute the Casimir force 
between a parabolic cylinder and a plate. In the limiting 
case when the radius of curvature at its tip vanishes, the 
parabolic cylinder becomes a semi-infinite plate (a knife's 
edge), and we can consider how the energy depends on 
the boundary condition it imposes and the angle it makes 
to the plane. 

The surface of a parabolic cylinder in Cartesian co- 
ordinates is described hy y = (a;^ — R^)/2R for all z, 



as shown in Fig. (TJ where R is the radius of curvature 
at the tip. In parabolic cylinder coordinates defined 
through X = /iA, y = (A^ — yU^)/2, z = z, the surface is 
simply fi = fj,o = \/R for — oo < A, z < oo. One advan- 
tage of the latter coordinate system is that the Helmholtz 
equation 

A2 + ^2 y gj^2 J dz^ ' V / 

which we consider for imaginary wavenumber k ~ in, 
admits separable solutions. Since sending A — — A and 
/i — >■ — /i returns us to the same point, we restrict our 
attention to /.t > while considering all values of A. Then 
jjL plays the role of the "radial" coordinate in scattering 
theory and we have regular and outgoing wave solutions 

^-(r) = e"''^D,{\)D.,^^{f,), (2) 

where D^{u) is the parabolic cylinder function, and 

A = A-y/ '2-\/k1 + and similarly for /i. Enforcing the 
reflection symmetry A — > — A and fi ^ —fi for the reg- 
ular solutions restricts the separation constant v to in- 
teger values. The corresponding outgoing solutions do 
not obey this restriction and thus can only be used away 
from /i = 0; as is typical for outgoing solutions, they are 
irregular at fi = 0. For imaginary wavenumber, the reg- 
ular (outgoing) solutions grow (decay) exponentially in 
/i and both i'^D^{iji) and Di,{\) are real. We can then 
express the free scalar Green's function as Q 

T^Et7^^^"(^<)*^^"'(^>)' 

(3) 
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where r< (r>) is the coordinate with the smaller (larger) 
value of /X. We will also use the Green's function in co- 
ordinates appropriate to scattering from a plane perpen- 
dicular to the y-axis, 



G(ri,r2,K) = / — e 



dk 



^ ik^{x2-xi)+iky \y2-y1 



.(4) 



where ky = i y/n^ + k^^+k^ . We can connect the 
parabolic and Cartesian Green's functions using the ex- 
pansion of a plane wave in regular parabolic solutions Q 



tan ■ 



/I 



-V'r(r), 



(5) 




FIG. 1: Parabolic cylinder/plane geometry. 



where tan (/) = This expression converges in regions 

where |tan^| < 1. A plane wave with |tan^| > 1 can 
instead be expanded in terms of solutions with negative 
integer values of [§] , and the Green's function can also 
be expressed in terms of these functions analogously to 
Eq. ([3]). Restricting to > is sufficient for our cal- 
culation, however, because we can already construct the 
Green's functions from these solutions alone; in the for- 
malism of Refs. [H, @, all possible quantum fluctuations 
are captured through the Green's function. Equating 
Eqs. ^ and ^ and then using ([S]), we obtain the expan- 
sion of the outgoing parabolic solution in plane waves. 



dhr 



I I tan I 



cos ■ 



-ikyy+ik^x 



, (6) 



which is valid for y < 0. 

The regular and outgoing waves provide two indepen- 
dent solutions to the second-order differential equation. 
We take a linear combination of these solutions to ob- 
tain the scattering solution $^(r) outside the parabolic 
cylinder. Fixing the coefficients by imposing Dirichlct 
boundary conditions at /i = /ig, we obtain 

$,(r) = Z?_„_i(/io)Vr(r) - t''D,{tf,o)rj'ir) , (7) 

while for Neumann boundary conditions we have 

= D'_,_,ii,o)^--ir) - i''+'Dliztlo)r.-^ir) . (8) 

These solutions to the Helmholtz equation can be used 
to compute the Casimir forces between a parabolic cylin- 
der and other simple objects, for example an infinite 
plate, as depicted in Fig. [TJ If both objects are perfect 
mirrors, translational symmetry along the z-axis enables 
us to decompose the electromagnetic field into two scalar 
fields, with Dirichlet and Neumann boundary conditions 
respectively. Each scalar field can then be treated in- 
dependently, with the sum of their contributions giving 



the full electromagnetic result. The quantization of each 
scalar field is achieved by integrating the exponentiated 
action over all configurations of the field [l^l • Constrain- 
ing the fields to obey the boundary conditions on each 
surface leads to an alternative description involving_fluc- 
tuating "charges" Ppia„c and p^^ 
Appropriate multipolcs of these charges are 



on the surfaces 



Q {kx k) 



dx dz dt e 



-ik^x — ik^z+Kt 



dX dz dt ( 



Ppia„c(a;,2:,t), 

ikz z+Kt 



y\^2mA J cylindc: 

/^o)/5cyi.„d„.(A, /io, z, t){\^ + pI). (9) 

The action can be decomposed in terms of these multi- 
poles as 5 = dn^^^^ [Spp + Sec + Scp + c.c], with 



Spp{K,k^) 

Scc{K,k:,) 
Scp{K,kz) 

00 

E 



—I 
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°° dk 
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r(k^r{Hx^Qp{k^), 
= -jEof (10) 



Here Spp corresponds to the action for the charges on 
the plane, with scattering amplitudes = ±1 for Neu- 
mann and Dirichlct modes respectively. The correspond- 
ing action for charges on the parabolic cylinder Sec 



can be related to its scattering amplitudes J-'^ 
Eqs. dl]) and ^ we obtain 



I; from 



■c 



^-.-i(Ao) 



(Dirichlct), 
(Neumann). 



(11) 



The position and orientation of the parabolic cylinder 
relative to the plane enter only through the translation 
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matrix lAyk^{d, which appears in the interaction term 
Sep- From Eq. (|6l), we obtain 



(tan 



2kyv\^/2^ cos2±i 



/ ^ikyd 



(12) 



where 6 is the angle of inchnation of the parabohc cylin- 
der and d is the distance from the focus of the parabola 
to the plane, as shown in Fig. [1] 

Integrating over these charge fluctuations gives the 
Casimir energy per unit length as 



£ 

hcL 



OO 1 

CLK 

2^ 

rOO 



2tt 



log det 1 



(13) 
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Numerical computations are performed by truncating the 
determinant at index i^„„„- For the numbers quoted be- 
low, we have computed for 1^,-^^^^ up to 200 and then ex- 
trapolated the result for — > oo. and in the figures 
we have generally used J^„ax = 100. We note that the 
integrals over k and fc^ can be expressed as a single in- 
tegral in polar coordinates, and for 6 = the integral 
is symmetric and the translation matrix elements vanish 
for v+h'' odd. Since the plane we are considering is a per- 
fect mirror, J-'^ is independent of k^ and we can further 
simplify the calculation for 6' = using the integral 



tan ■ 



2n 



' ky cos2 I 



2^k_2„-i(2dVK2 + fc2), 



(14) 

where k£{u) = —fj — -^C/(— |, 0, 2u) is the Bateman k- 

n > 
function [llj, which is zero if £ is a negative even integer. 

Here U{a,b,u) is the confluent hypergeometric function 

of the second kind. 

As a first demonstration, we report on the dependence 

of the energy on the separation H = d — R/2 ior 6 = 0. 

At small separations [H / R <C 1) the PFA, given by 



hcL 



720 



dx 



-OO [H + x^/{2R)Y 



960V2 



(15) 

should be valid. The numerical results in Fig. [5] confirm 
this expectation with a ratio of actual to PFA energy 
of 0.9961 at H/R = 0.25 (for i? = 1). We note that 
since the main contribution to PFA is from the proximal 
parts of the two surfaces, the PFA result in Eq. (fT5|) also 
applies to a circular cylinder with the same radius R. 

A more interesting limit is obtained when RjH — > 0, 
corresponding to a semi-infinite plate. Then the PFA 
result is zero, as are results based on perturbative ap- 
proximation for the dilute limit 12] ■ The scattering am- 



plitudes in Eq. (jll[) simplify and can be combined to- 



gether as T\ 



■c 



'2/tt, where even v corresponds to 
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FIG. 2: The energy per unit length times H'^ , EH^/(hcL), 
plotted versus H/R for 6 = and i? = 1 on a log-linear scale. 
The dashed line gives the R = limit and the solid curve 
gives the PFA result. 



Dirichlet and odd u corresponds to Neumann. Using this 
result, our expression for the energy for R = and 9 = 
simplifies to 



£ 

hcL 



^°°^logdet(l,,. - (-l)''k_,_,,_i(2gi?)) 

(16) 



where C±^ = 0.0067415 is obtained by numerical inte- 
gration. This geometry was studied using the world-line 
method for a scalar field with Dirichlet boundary con- 
ditions in Ref. [13|- The world- line approach requires a 



large-scale numerical computation, and it is not known 
how to extend this method to Neumann boundary con- 
ditions (or any case other than a scalar with Dirichlet 
boundary conditions). In our calculation, the Dirichlet 
component of the electromagnetic field makes a contribu- 
tion = 0.0060485 to our result, in reasonabl e ag ree- 
mcnt with the value of Cf = 0.00600(2) in Ref. [IJ. 
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FIG. 3: The coefficient c{6) as a function of angle for R = Q. 
The exact result at = 7r/2 is marked with a cross. Inset; 
Dirichlet (circles) and Neumann (squares) contributions to 
the full electromagnetic result. 
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Reference [l3| also considers a tilted semi-infinite plate, 
which corresponds to the i? — >■ limit of our formula for 
general 9. From dimensional analysis, the Casimir energy 
at i? = again takes the now 0-dependent form 



£ 

hcL 



(17) 



where H ~ d for R = 0. Following Ref. [13[, we 
plot c{e) = cos{0)C{9) in Fig. El A particularly in- 
teresting limit is 6 — >■ 7r/2, when the two plates are 
parallel. In this case, the leading contribution to the 
Casimir energy should be proportional to the area of 
the half-plane according to the parallel plate formula, 
E^^/{ficA) = -C||/i7^ with C|| = 7r^/720, plus a subleading 
correction due to the edge. Multiplying by cos 9 removes 
the divergence in C{9) as — > 7r/2. As in Ref. we 
assume c{9 — )- 7r/2) = c\\/2 + {9 — 7r/2) Cedgo, although we 
cannot rule out the possibility of additional non-analytic 
forms, such as logarithmic or other singularities. With 
this assumption, we can estimate the edge correction 
Codgo = 0.0009 from the data in Fig. [31 From the inset in 
Fig. 131 we estimate the Dirichlet and Neumann contribu- 



tions to this result to be c 



-0.0025 (in agreement 



N 



0.0034 



with [13l | within our error estimates) and c 
respectively. Because higher partial waves become more 
important as — >■ 7r/2, reflecting the divergence in C{9) 
in this limit, we have used larger values of ly^^,, for 9 near 
7r/2. 

It is straightforward to extend these results to nonzero 
temperature T. We simply replace the integral |^ 
by the sum -^X^J^o' "^^^r Matsubara frequencies k„ = 
2TmT/{hc), where the prime indicates that the n — 
mode is counted with a weig ht of 1/2 0. In the limit 
of infinite temperature, only the n = mode contributes 
and we obtain for R = the energy £/L — ~TCt=oo/ H , 
with Ct=oo = 0.0472. The Dirichlet contribution to our 
result is C|Lqq = 0.0394, again in agreement with [is} . 

Employing the scattering formalism, we can also cal- 
culate the Casimir energy for the case where another ob- 
ject whose scattering amplitudes are available, such as 
an ordinary cylinder or a second parabolic cylinder, is 
positioned outside the parabolic cylinder. Centering the 
other object at the origin and letting the parabolic cylin- 
der open downward, with its focus displaced to y = —d, 
we obtain the necessary translation matrix elements by 
writing Eq. ([6]) for f, where x — x, y ~ ~y — d, z = z, 
and then expanding the plane wave on the right-hand 
side in the basis appropriate to the other object. Again 
we can allow the parabolic cylinder to tilt by replacing cf) 
hy cl)+9 in this expression. These results can be extended 
to multiple objects, as in Ref. [lj|. Another interesting 
possibility would be to apply the interior Casimir for- 
malism of Ref. [isl l an object inside a parabolic cylinder, 
potentially extending the results of Ref. [3, 12 1. 

The reduction of the parabolic cylinder to a semi- 
infinite plate enables us to consider a variety of edge ge- 



ometries. A thin metal disk perpendicular to a nearby 
metal surface would experience a Casimir force described 
by an extension of Eq. (fT6|) . Figure [21 shows that the 
PFA breaks down for a thin plate perpendicular to a 
plane; the PFA approximation to the energy vanishes 
as the thickness goes to zero, while the correct result 
instead has a different power law dependence on the 
separation. Based on the full result for perpendicular 
planes, however, we can formulate an "edge PFA" that 
yields the energy by integrating d£/dL from Eq. (flB]) 
along the edge of the disk. Letting r be the disk radius, 
in this approximation we have ^Epta = —ficC± J_,^{H + 



x^) '^dx Jllll^ 



-hcC±n.yr/{2H^), which is 
valid if the thickness of the disk is small compared to its 
separation from the plane. (For comparison, note that 
the ordinary PFA for a metal sphere of radius r and a 
plate is proportional to r/H'^.) 

A disk may be more experimentally tractable than a 
plane, since its edge does not need to be maintained par- 
allel to the plate. One possibility would be a metal film, 
evaporated onto a substrate that either has low permit- 
tivity or can be etched away beneath the edge of the de- 
posited film. Micromechanical torsion oscillators, which 
have already been used for Casimir experiments [l8| . 
seem readily adaptable for testing Eq. pT]) . Because the 
overall strength of the Casimir effect is weaker for a disk 
than for a sphere, observing Casimir forces in this geome- 
try will require greater sensitivities or shorter separation 
distances than the sphere-plane case. As the separation 
gets smaller, however, the dominant contributions arise 
from higher-frequency fluctuations, and deviations from 
the perfect conductor limit can become important. While 
the effects of finite conductivity could be captured by an 
extension of our method, the calculation becomes signif- 
icantly more difficult in this case because the matrix of 
scattering amplitudes is no longer diagonal. 

To estimate the range of important fiuctuation fre- 
quencies, we consider i? <C -ff and 9^0. In this 
case, the integrand in Eq. ([T5| is strongly peaked around 
q « 0.3/H. As a result, by including only values of q up 
to 2/H, we still capture 95% of the full result (and by 
going up to 3/iJ we include 99%). This truncation corre- 
sponds to a minimum fluctuation wavelength An,i„ = ttH. 
For the perfect conductor approximation to hold, A„,i„ 
must be large compared to the metal's plasma wavelength 
Xp, so that these fluctuations are weU described by as- 
suming perfect reflectivity. We also need the thickness 
of the disk to be small enough compared to H that the 
deviation from the proximity force calculation is evident 
(see Fig. [2]), but large enough compared to the metal's 
skin depth 6 that the perfect conductor approximation 
is valid. For a typical metal film, Xp « 130 nm and 
S Ri 25 nm at the relevant wavelengths. For a disk of 
radius r = 100 ^m, the present experimental frontier of 
0.1 pN sensitivity corresponds to a separation distance 
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H « 350 nm, which then falls within the expected range 
of validity of our calculation according to these criteria. 
The force could also be enhanced by connecting several 
identical but well-separated disks. In that case, the same 
force could be measured at a larger separation distance, 
where our calculation is more accurate. 
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